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The  effect  of  shear  stress  on  the  magnitude  and  directional  variation 
of  the  dynamic  bulk  modulus  of  sand  was  investigated  by  means  of  a program 
of  laboratory  measurements  of  P-wave  and  S-wave  velocities  under  varied 
stress  conditions.  The  sand  used  was  air-dry  Reid-Bedford  Model  sand  and 
was  tested  at  relative  densities  of  about  25%  and  80)C . Test  specimens  were 
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of  Octahedral  normal  stress  of  5,  10,  and  20  psi,  and  with  ratips  of 

lyWid  1/3.  ^Measurement^  of  P-wave  and  S-wave  velocities  were  made  in  three  or 
four  directions  with  Accelerometers  buried  in  the  test  specimens. 

— > The  test  results  indicate  that  compression  and  shear  wave  velocities 

vary  with  direction  in  an  isotropic  stress  field  due  to  inherent  anisotropy 
of  the  sand  structure , hut  that  this  phenomenon  has  only  a modest  ( 0-20% ) 
effect  in  producing  anisotropy  in  the  bulk  modulus.  The  work  tends  to  confirm, 
within  +10){,  that  the  shear  wave  velocity  depends  only  on  the  level  of  octahe- 
dral normal  effective  stress  and  not  on  stress  anisotropy.  In  contrast  with 
the  shear  wave  velocity,  the  compression  wave  velocity  increases  significantly 
in  the  direction  of  the  major  principal  stress  and  depends  primarily  on  the 
effective  stress  on  the  plane  normal  to  wave  travel.  Values  of  bulk  modulus, 
computed  from  the  wave  velocities  under  the  assumption  of  elastic  isotropy, 
show  a variation  with  direction,  which  contradicts  the  assumption  of  isotropy. 

Appendix  A contains  photographs  of  typical  oscilloscope  waveforms  from 
each  replicate  series  of  five  impacts,  and  Appendix  B presents  tables  of  cor- 
rected V and  V , the  push  wave  and  shear  wave  velocities,  respectively, 
from  each^test. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 


verted to  metric  (Si)  units  as  follows: 


Multiply 

cubic  feet 
degrees  (angle) 


feet  per  second 
inches 

pounds  (force)  per 
square  inch 

pounds  (mass) 

pounds  (mass)  per 
cubic  foot 


0.02831685 

0.171+5329 

0.3048 

0.3048 

25.4 

6894.757 

0.4535924 

16.01846 


To  Obtain 

cubic  metres 

radians 

metres 

metres  per  second 

millimetres 

pascals 


kilograms 

kilograms  per  cubic 
metre 
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EFFECT  OF  SHEAR  STRESS  ON  DYNAMIC  BULK  MODULUS  OF  SAND 


1.  INTRODUCTION 

1.1  Purpose 

The  hulk  modulus,  B,  represents  one  of  the  most  important  soil 
properties  that  determine  the  possible  liquefaction  behavior  of  a cohe- 
sionless soil.  The  greater  the  value  of  B,  the  lower  the  magnitude  of 
pore  pressure  response  to  cyclic  loading  and  the  less  likely  is  lique- 
faction. 

In  accord  with  the  theory  of  elasticity  in  a homogeneous,  iso- 
tropic medium,  most  if  not  all  investigators  have  assumed  B is  an 
isotropic  property  independent  of  the  direction  of  the  travel  of  the 
disturbing  stress  waves  that  might  induce  liquefaction.  The  new  theo- 
retical analysis  proposals  by  F.  E.  Richart,  Jr.,  and  co-workers  at  the 
University  of  Michigan  have  caused  new  questions  to  arise  about  the 
validity  of  the  assumption  of  isotropic  B.  Accordingly,  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES),  CE,  initiated  research  to 
make  at  least  a preliminary  assessment  of  the  degree  of  validity  of  the 
isotropic  B assumption. 

In  an  isotropic,  homogeneous,  elastic  medium: 

B=‘>(Vp-IVs)  (1) 

Thus,  we  need  only  determine  the  push  (compression)  wave  and  shear  wave 

velocities,  V and  V , respectively,  together  with  the  mass  density 
P ® 

of  the  sand,  p , to  determine  its  bulk  modulus  when  subjected  to  stress 
wave  loadings.  We  know  that  sands  behave  as  particulate  systems  and 
that  many  of  their  properties  depend  on  the  level  of  ambient  effective 
stress.  Previous  research  (Hardin  and  Black,  1966)  showed  that  the 
shear  wave  velocity  depended  only,  within  about  +5%,  on  the  magnitude 
of  octahedral  effective  stress  in  a soil  element  and  not  on  the  ratio 
of  the  principal  stresses  and  therefore  the  shear  stress  magnitudes. 
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This  investigation  had  the  primary  purpose  of  determining  if  B also 
depended  only  on  the  magnitude  of  octahedral  stress  and  to  check  its 
independence  from  the  direction  of  wave  propagation  when  in  an 
anisotropic  stress  field. 


1.2  Scope 


The  scope  of  this  work  includes  a single  sand,  tested  air-dry, 
following  pluvial  deposition  to  two  relative  densities,  about  25 % and 
80%,  corresponding  to  dry  density,  y^,  of  93  and  103  pcf.*  Four  sensi- 
tive accelerometers,  buried  in  the  sand  during  the  casting  of  each  of 
ten  separate  specimens,  permitted  determination  of  the  V and  Vg 
velocities  in  three  or  four  directions  within  a U-ft-diam,  U-ft-high 
sand  specimen  in  a test  chamber.  At  each  relative  density  condition, 
we  performed  a test  with  K = o^/o^  = 1 and  K = 1/3  to  evaluate  the 
importance  of  nonequal  principal  stresses  on  the  behavior  of  V , Vg , 
and  B for  the  different  directions  of  wave  travel.  For  each  test  and 
for  each  direction  of  wave  travel,  we  made  out  measurements  at  three 
levels  of  constant  octahedral  stress  a ^ — 5»  10,  and  20  psi.  Further- 
more, we  performed  either  two  or  three  replicate  tests  at  each  of  the 
relative  density  and  K conditions. 


In  sum,  we  made  over  five  hundred  separate  determinations  of  each 
and  Vg.  Our  summarized  results  show  a minor,  perhaps  negligible, 
dependence  of  B on  wave  propagation  direction  under  isotropic  (K  = l) 
stress  conditions,  but  a strong  dependence  of  B on  direction  when  sub- 
jected to  an  anisotropic  (K  = 1/3)  stress  condition. 


2 . EQUIPMENT 


2.1  Test  Chamber 

The  key  piece  of  equipment  that  made  this  research  possible  was 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (SI)  units  is  given  on  page  5. 


the  4-ft-diam,  4-ft-high,  triaxial  chamber  specimen  within  which  we 
placed  accelerometers  during  the  pluvial  deposition  used  to  cast  the 
specimens.  This  University  of  Florida  calibration  chamber  permits 
independent  control  of  the  radial  and  vertical  boundary  stresses 
on  the  specimen.  The  pluvial  deposition  permits  the  casting  and  test- 
ing of  approximately  duplicate  specimens.  The  sand  falls  from  a travel- 
ing hopper,  through  shutter  plates,  through  two  layers  of  coarse  screen, 
and  then  into  the  chamber.  By  varying  the  size  of  the  holes  in  the 
shutter  plate,  we  can  achieve  different  relative  densities.  Figure  1 
presents  a schematic  cross  section  through  this  test  chamber.  Please 
refer  to  other  references  ( Schmertmann , 1978:  Veismanis,  1974:  Laier 
et  al.,  1975)  for  more  detailed  descriptions  of  the  design  and  opera- 
tion of  this  chamber. 

The  work  of  Caillemer  (1975)  and  others  has  shown  that  a given 

shutter  plate  hole  size  produces  repeatable  average  chamber  densities 

within  +1  pcf , and  that  dry  unit  weight  variations,  both  vertically 

3 

and  radially  within  the  50-ft  test  specimen,  vary  randomly  within 
+1 . 5 pcf  from  the  whole  specimen  average.  Caillemer  tested  with 
samples  having  a volume  of  492  cm  , or  0.035%  of  the  total  specimen,  and 
took  20  samples  per  specimen. 

2.2  Sand 

For  all  tests  in  this  research  we  used  only  one  sand,  the  Reid- 
Bedford  Model  sand.  We  obtained  this  sand  from  WES  in  1976.  Multiple 
use  of  the  sand  through  our  pluvial  deposition  equipment  and  subsequent 
crop-type  dryer  equipment  in  those  cases  in  previous  research  where  we 
wetted  the  sand  has  resulted  in  our  losing  most  of  the  minus  200  sieve 
fines.  Figure  2 shows  the  results  of  several  recent  sieve  analyses  on 
the  sand  we  used. 

2.3  Accelerometers  and  Oscilloscope 

In  each  of  the  chamber  test  samples,  we  buried  four  acceler- 
ometers, all  Piezotronics  model  No.  308A02,  with  a sensitivity  of 
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100  mv/g.  Fach  of  these  accelerometers  had  its  own  battery-operated 
signal  conditioning  unit  supplied  by  the  manufacturer.  In  general,  the 
accelerometers  appeared  to  work  well  for  all  tests.  We  did  have  some 
trouble  with  the  conditioning  units,  which  in  the  final  few  tests 
required  the  switching  of  accelerometers  A1  and  A2  (numbering  discussed 
subsequently)  into  the  conditioning  units  for  A3  and  A*4.  While  this 
affected  the  calibration  for  acceleration  magnitude,  it  did  not  influence 
the  more  important  time  intervals  between  wave  arrivals. 

By  the  use  of  an  appropriate  triggering  system,  we  displayed  the 
response  from  each  pair  of  accelerometers  on  a dual,  chopped-beam , 
storage  oscilloscope.  We  used  an  HP  l8lA  scope  with  a 180TA  dual 
channel  vertical  amplifier  and  a 1821A  time  base  generator.  We  used 
only  the  0.2,  0.5,  and  1.0  msec/division  (ms/d)  sweep  rates  for  the 
key  determinations  of  the  differences  in  the  times  of  first  arrival  for 
the  push  and  shear  waves  at  each  accelerometer.  After  the  research,  we 
calibrated  these  sweep  rates  and  corrected  the  velocity  data  for  the 
small  (about  0.5  percent)  errors  involved. 


Figure  3,  photo  (a),  shows  the  transducers  we  used.  Photo  (b) 
shows  the  oscilloscope.  Photo  (c)  shows  two  of  the  accelerometers, 
Nos.  2 and  3,  being  carefully  laid  at  their  positions  for  Tests  6 
through  13  during  the  filling  of  the  sand  in  the  chamber.  Photo  (d) 
shows  the  operators  monitoring  the  oscilloscope  during  the  performance 
of  the  test. 


2.U  Impact  Equipment 


For  Tests  3 and  4,  we  used  a hand-held  and  swung  ball  peen  hammer 
for  our  wave-producing  impact  system.  It  soon  became  evident  that  we 
could  improve  operations  significantly  by  introducing  a mechanical  im- 
pact system  to  make  our  impacts  and  resulting  waveforms  less  operator- 
dependent  and  therefore  more  reproducible.  The  photos  in  Figure  U show 
the  mechanical  system  we  used  for  Tests  7 through  13.  (Test  5 was 
aborted  when  the  need  for  a better  impact  system  became  evident.) 


The  mechanical  impact  system  consisted  of  either  dropping  or 
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swinging  a steel  ball  weighing  176  grams.  When  we  wanted  to  produce 
waves  relatively  rich  in  compression  energy,  we  dropped  the  ball  from 

2 to  6 in.  onto  the  center  of  a steel  plate  capping  the  vertical  steel 
rod  extending  from  aluminum  plates  buried  at  two  positions,  each 

3 cm  below  the  top  of  the  sand  specimen.  When  we  wanted  to  produce 
waves  relatively  rich  in  shear  energy,  we  swung  the  weight  as  a pendulum 
and  struck  a horizontal  bar  attached  to  a horizontally  braced  (with 
plastic  pipe)  vertical  steel  rod  from  the  buried  plate.  Both  types  of 
impact  involved  kinetic  energies  at  impacts  of  about  1.0  kg-cm,  +50$. 

3.  METHODS 


3.1  Placement  of  Accelerometers 

As  noted  above,  and  also  illustrated  in  Figure  3(c),  we  placed 
the  accelerometers  in  the  chamber  during  appropriate  interruptions  in 
the  pluvial  filling  process.  We  placed  all  accelerometers  in  a single, 
vertical  plane  passing  through  the  chamber  axis.  When  the  sand  filling 
reached  an  appropriate  depth,  we  put  each  accelerometer  in  its  planned 
location  with  an  estimated  precision  of  better  than  +1/U  cm  in  both 
depth  and  radial  position.  After  some  initial  experimentation  with 
alternate  orientations  of  the  axes  of  the  accelerometers,  I decided  to 
place  all  accelerometers  with  their  axes  horizontal  and  perpendicular 
to  plane  in  which  we  placed  all  the  accelerometers.  All  the  accelerom- 
eters had  their  wires  coming  out  of  the  same  side,  giving  all  the  same 
polarity. 

Figure  5 shows  my  initial  choice  for  the  position  of  the  four  ac- 
celerometers along  a vertical  plane  cutting  through  the  chamber  diameter 
and  through  the  center  of  the  two  impact  access  holes  in  the  top  of  the 
chamber.  The  accelerometer  positions  shown  in  Figure  5 are  the  ones  we 
used  for  Tests  3 and  h.  These  positions  had  certain  advantages,  such 
as  providing  three  wave  speed  measurement  directions  at  the  same  average 
magnitude  of  octahedral  stress.  However,  we  found  serious  interference 
at  accelerometers  A3  and  Ah  due  to  their  position  near  the  chamber 
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boundaries  and  the  consequent  first  arrival  of  the  waves  traveling 
through  the  boundary  steel  and  water.  To  help  correct  this,  I changed 
to  a new  pattern  for  the  accelerometers,  shown  in  Figure  6,  which  we 
used  for  Tests  6 through  13.  This  pattern  moved  the  outermost  acceler- 
ometers farther  away  from  the  walls  of  the  chamber  and  provided  for 
velocity  measurements  in  four  directions  with  the  horizontal,  90°  (A2- 
AU),  60°  (A3-AU),  30°  (A1-A2),  a 0°  (A3-A2).  The  Figure  6 pattern 
proved  reasonably  successful , although  we  did  not  completely  eliminate 
chamber  boundary  wave  interference  effects  at  the  furthermost  accelerom- 
eter in  each  paired  direction. 

We  placed  each  accelerometer  in  position  with  a minimum  of  bedding 
disturbance  to  the  surrounding,  pluvially  placed  sand,  with  the  wires 
going  directly  to  the  side  of  the  chamber  away  from  the  plane  of  the 
velocity  testing.  These  wires  eventually  came  out  through  the  central 
hole  in  the  chamber,  but  we  kept  them  as  far  as  possible  from  the  plane 
of  the  testing.  We  placed  each  accelerometer  with  its  protective  rubber 
cap  as  recommended  by  the  manufacturer.  Figure  3(e)  shows  a photo  of 
one  of  these  accelerometers  with  the  cap  removed. 

3.2  Applying  the  Chamber  Pressure 

For  each  test,  we  controlled  the  chamber  pressure  by  means  of  two 
pressure  regulators,  one  controlling  vertical  pressure  and  one  control- 
ling radial  pressure.  The  entire  bottom  of  the  specimen  rested  on  a 
hydraulic  load  cell,  constructed  as  a water-filled,  neoprene-covered, 
disc-shaped  recess  in  the  top  of  the  vertical  piston.  This  piston 
acted  at  the  bottom  of  the  chamber  and  specimen.  Attached  to  this  load 
cell  on  top  of  the  piston  was  a Heise  Bourdon  gauge  with  0.2-psi  divi- 
sions and  a readability  of  +0.05  psi,  and  connected  to  the  water  a- 
round  the  circumference  of  the  specimen  was  a similar  Heise  Bourdon 
gauge  to  monitor  the  water  pressure  in  this  circumf erential  water. 

During  the  application  of  the  sample  confining  pressure,  the 
operator  used  air  pressure  regulators  and  the  Heise  gauges  to  increase 
both  radial  and  vertical  pressures  in  small  increments,  according  to  a 
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set  constant-K  path  that  depended  on  whether  the  test  called  for  K = 1 
or  K = 1/3  conditions.  Except  for  an  initial  adjustment,  he  maintained 
the  proper  K loading  path  until  reaching  the  specified  test  levels  of 
5.0-,  10. 0-,  and  20.0-psi  average  octahedral  stress  midheight  between 
accelerometers  A2-AH  and  A3-AU. 

We  applied  chamber  pressure  only  in  the  direction  of  increasing 
stress.  None  of  the  samples  had  ever  been  subjected  to  higher  levels 
of  effective  stress  than  they  had  at  the  time  of  testing. 

3.3  Production  of  the  Impact  Waves 

For  this  purpose  we  used  the  impact  equipment  described  in  section 
2.h.  With  the  mechanical  system,  the  ball  impacted  against  steel,  and 
the  resulting  stress  wave  traveled  down  the  vertical,  standard  pipe  rod 
of  1/2-in.  diameter.  This  steel  rod  was,  in  turn,  rigidly  attached  to 
a horizontal  aluminum  plate  or  disc  of  3-in.  diameter  and  l/l+-in.  thick- 
ness. We  buried  this  plate  with  its  bottom  approximately  3 cm  below  the 
top  of  the  specimen.  We  had  two  such  plates,  one  below  the  center  ac- 
cess hole  and  one  below  the  side  access  hole.  The  bottom  of  each  plate 
had  rough  sandpaper  glued  to  it  to  assure  good  frictional  contact  be- 
tween soil  and  plate.  The  vertical  rods  coming  up  from  the  plate  were 
also  braced  with  concentric  plastic  pipe.  This  pipe  had  the  dual  pur- 
pose of  electrically  isolating  the  vertical  rod  from  the  chamber  and 
thus  preventing  early  triggering  of  the  scope  impact  sweep  and  also  of 
providing  a reaction  against  the  rocking  component  of  the  impact  gener- 
ated by  the  horizontal,  pendulum  ball  impact. 

With  the  system  we  used,  the  operator  would  clear  the  storage 
scope,  set  it  for  the  next  sweep,  and  then  call  for  the  impact  device 
operator  to  either  drop  or  allow  the  ball  to  swing  and  make  the  impact. 
The  contact  of  the  ball  with  the  vertical  rod  system  would  then  complete 
the  triggering  circuit  and  start  the  sweep  of  the  scope  beams  at  approxi- 
mately the  same  instant  as  the  stress  wave  started  in  the  rods,  through 
the  plate,  and  then  through  the  sand  to  the  accelerometers.  We  used  a 
separate  5-v  system  to  actuate  the  triggering  for  all  tests. 
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For  Tests  9 through  13,  we  also  measured  V along  the  horizontal 

P 

direction  between  accelerometers  A3-A2  (see  Figure  6).  To  do  this,  we 
used  the  mechanical  ball  system,  but  hand-held  the  ball  and  struck  it 
against  the  steel  outside  of  the  tank  at  the  point  where  the  projection 
from  A2-A3  intersects  the  side  of  the  tank.  Despite  attempts,  we  could 
not  interpret  the  resulting  waveforms  for  V . 

3.^  Interpreting  and  Recording  the  Impact  Waveforms 

Following  the  impact,  the  oscilloscope  displayed  the  waveform 
record  as  sensed  by  each  pair  of  accelerometers,  displaying  this  by  use 
of  the  chopped  mode.  The  operator  at  the  scope  (usually  the  writer) 
would  call  for  a succession  of  impacts  until,  by  manipulation  of  the 
gain  and  sweep  rate  controls,  he  obtained  a reproducible  and  informative 
waveform.  First  he  called  for  vertical  (dropped  ball)  impacts  to  get 
the  first  arrival  times  of  the  push  (compression)  waves  at  any  two  of 
the  matched  accelerometer  locations.  The  actual  arrival  times  at  each 
were  only  of  concern  with  respect  to  the  reasonableness  of  the  waveform 
interpretation  for  first  push  wave  arrival.  Of  greatest  concern  was  an 
accurate  measure  of  the  time  interval  between  the  first  arrival  at  the 
two  accelerometers.  To  obtain  this  difference,  the  operator  made  an 
estimate,  directly  from  the  stored  scope  waveform  record,  of  the  time  of 
first  arrival  to  +0.05  major  time  divisions  on  the  face  of  the  scope. 

To  achieve  better  accuracy  in  the  final  wave  velocity  computation,  we 
repeated  the  above-mentioned  impacts  until  we  obtained  five  waveforms 
that  we  could  interpret  for  time  differences  in  the  first  arrival.  These 
five  were  called  out  to  a recorder  (usually  a graduate  student).  We 
then  averaged  these  five  for  the  time  difference  used  to  calculate 

V (and  V ) . 

P s 

The  scope  operator  then  repeated  the  above-mentioned  procedure 
with  the  swinging  ball  for  an  interpretation  of  the  first  arrival  times 
of  the  shear  wave.  This  usually  involved  a more  difficult  interpretation 
because  the  arrival  of  the  shear  wave  was  always  superposed  on  the  now 
already  present  push  wave.  Unfortunately,  with  our  impact  system  we 
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were  not  able  to  obtain  a reversal  in  the  polarity  of  the  shear  wave 
without  a simultaneous  reversal  in  the  polarity  in  the  accompanying 
push  waves.  Thus,  reversal  did  not  aid  interpretation.  Our  oscillo- 
scope did  not  have  an  enhancement  capability.  However,  we  doubt  that 
this  would  have  improved  the  ease  of  interpretation  because  we  had 
virtually  no  random  background  noise  with  which  to  contend. 

3.5  Photographing  the  Impact  Waveforms 

During  each  series  of  the  five  scope  displays  of  waveforms  for 
which  we  interpreted  first  arrival  times  at  each  accelerometer,  the 
operator  chose  a typical  set  of  waveforms  and  photographed  them,  using 
a suitable  Polaroid  camera.  The  recorder  noted  the  particular  set  of 
readings  for  which  a photograph  was  taken.  The  operator  then  im- 
mediately noted  on  the  photograph  its  sequence  number,  and  such  other 
information  as  test  number,  test  time,  test  pressures,  and  scope  gain 
and  sweep  settings,  and  most  important — his  assumption  for  the  push 
and/or  shear  wave  arrival  times  for  that  waveform.  Except  for  a few 
readings  in  Test  3,  all  the  series  of  five  waveforms  have  at  least  one 
photograph  documenting  a typical  waveform  and  our  interpretation. 

Appendix  A presents  annotated  copies  of  all  the  waveform  photo- 
graphs we  took  to  document  this  research.  At  the  time  of  taking  the 
photograph  each  was  given  the  next  successive  acquisition  number.  The 
photographs  are  presented  in  order  of  acquisition  number,  given  in  the 
upper  right  corner  of  each. 

3.6  Measured  Average  Dry  Unit  Weight 

We  weighed  all  the  sand  removed  from  the  chamber  after  completion 
of  each  test.  We  knew  the  volume  of  the  specimen  at  the  beginning  of 
each  test  because  we  knew  the  initial  volume  of  the  chamber  (50.3  ft  ) 
and  the  volume  changes  measured  during  the  application  of  pressures 
before  each  test.  Thus,  we  could  easily  compute  the  average  air-dry 
unit  weight  of  the  sand  for  each  test.  Considering  the  accuracy  of 
our  spring-scale  weighing  system  (readability  lb)  and  the  number  of 


weighings  to  complete  weighing  all  the  sand  (about  5),  we  believe  we 
obtained  the  total  weight  of  the  sand  within  +25  lb.  In  terms  of  unit 
weight,  this  represents  about  +0.5  pcf.  In  two  tests  (3  and  8),  we  had 
a greater  error  because  some  of  the  sand  got  wet  due  to  membrane  leak- 
age prior  to  completion  of  the  wave  velocity  testing. 

3.7  Time  Sequence 

Table  1 presents  a chronological  listing  of  the  major  operations 
involved  in  one  of  our  tests,  with  the  purpose  of  giving  the  reader  a 
better  concept  of  the  sequence  of  events  in  a test  and  the  various 
typical  times  required. 

From  Table  1 the  reader  can  see  that  we  could  perform  between  one 
and  two  tests  per  week  and  that  each  test  required  approximately  45  man- 
hours of  work  to  complete,  independent  of  any  subsequent  plotting  and 
analyses  of  data. 


ANALYSES  OF  DATA 


4.1  Correction  of  Raw  Data 


As  explained  in  the  previous  sections,  the  basic  data  from  this 

research  consist  of  differences  in  the  arrival  times  for  the  push  and 

shear  waves,  which  permit  the  calculation  of  V and  V . Some  of  the 

p s 

associated  measurement  errors  are  random,  and  some  systematic.  In  turn, 
some  of  these  systematic  errors  can  be  estimated  and  corrected.  We  have 
included  such  corrections.  The  following  explains  what  we  did. 

4.1.1  Different  Levels  of  Octahedral  Effective  Stress  Between 
Ac c elerometer s . We  do  not  have  uniform  effective  stress  conditions  in 
a test  chamber  sample  because  of  the  effect  of  the  unit  weight  of  the 
sample  itself.  The  higher  the  point  in  the  sample,  the  lower  both  the 
horizontal  and  vertical  effective  stresses.  Reference  to  Figure  6 will 
show  that  the  average  effective  stresses  between  accelerometers  A1-A2 
and  A3-A2  must  be  less  than  between  A2-A4  and  A3-A4.  Because  both  push 
and  shea  wave  velocities  increase  with  effective  stress,  we  corrected 
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the  velocity  data  between  A1-A2  and  A3-A2  (angles  30°  and  0°)  to  make 
them  conform  to  the  same  5.0-,  10. 0-,  or  20.0-psi  octahedral  stress 
magnitude  between  accelerometers  A2-A4  and  A3-A4  (the  90°  and  60° 
directions) . 

The  formula  by  Hardin  and  Drnevich  (1972)  indicates  the  shear  wave 

velocity  will  vary  with  the  1/4  power  of  the  octahedral  stress.  A 

recent  unpublished  Japanese  reference  indicates  that  the  shear  wave 

velocity  varies  with  a little  less  than  the  1/5  power.  Plotting  some 

of  our  results  on  log-log  paper  indicates  that  both  V and  V vary  with 

P s 

about  the  1/5  power.  We  used  the  1/5  power  as  our  basis  for  increasing 
the  wave  velocities  measured  between  accelerometers  A1-A2  and  A3-A2  to 
have  these  velocities  conform  to  the  stated  octahedral  stress  level  of 
either  5»  10,  or  20  psi.  Doing  this  introduced  a maximum  correction  of 
+3 # at  the  5-psi  level,  reducing  to  a maximum  of  +0.07#  at  the  20-psi 
level . 


4.1.2  Compression  of  Sample  Under  Pressure.  With  the  application 
of  pressure  to  the  various  octahedral  stress  levels  for  the  test,  the 
sample  usually  changed  both  vertical  and  radial  dimensions.  Because  we 
measured  volume  changes  in  both  directions,  we  could  compute  the  average 
vertical  and  radial  strains.  By  assuming  these  strains  occur  uniformly 
throughout  the  sample,  we  corrected  the  original  as-placed  distances 
between  accelerometers  to  slightly  modified  distances  after  the  sample 
compressed  to  reach  each  successive  octahedral  stress  level.  Because 
the  distances  shorten,  this  correction  reduces  the  velocities  calculated 
by  use  of  the  as-placed  distances.  We  also  applied  this  correction, 
which  had  a typical  value  of  about  -1/2#  and  a maximum  of  -1#.  The 
angles  changed  a negligible  amount,  and  we  disregarded  this  effect. 


4.1.3  Unit  Weight  Variations.  For  the  purposes  of  comparing  the 

V and  V determinations  within  each  group  of  low  and  high  relative 
P s 

density  tests,  we  also  corrected  for  the  probable  differences  in  dry 
density  between  each  individual  test  and  the  average  for  the  entire 
group.  We  based  this  correction  on  the  measured  changes  in  and  Vg 
over  the  much  larger  increment  in  unit  weight  between  the  high  and  low 
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density  tests.  Over  this  increment  we  found  an  average  increase  in 
velocity  of  1.0#  for  each  1.0-pcf  increase.  We  then  used  this  as  a 
"basis  to  further  correct  the  wave  velocity  results  to  have  each  group 
conform  to  the  average  dry  unit  weight  for  that  group  (as  determined  by 
the  weighing  method  described  in  section  3.6). 

4.1.4  Scope  Sweep  Calibration.  Several  days  after  completion  of 
Test  13,  we  took  the  oscilloscope  used  in  this  research  to  the  Univer- 
sity of  Florida  Instrument  Standards  Laboratory  and  calibrated  the  sweep 
rates  with  a Tektronix  type-l8U  signal  generator.  We  only  used  the 
0.2,  0.5,  and  1.0  ms/d  sweep  rates  in  this  research.  In  all  cases  we 
found  the  sweep  rate  in  error  by  less  than  1%. 


As  a result  of  this  calibration,  we  used  the  following  factors 
by  which  to  multiply  and  Vg  to  correct  for  instrument  sweep  rate 
error:  0.9 94  for  the  0.2  scale,  and  0.996  for  both  the  0.5  and  1.0  ms/d 
scales. 


4.2  Graphs  and  Tables  of  Corrected  Velocities 

Figures  7 through  12  present  polar  coordinate,  graphical  sum- 
maries of  the  above-corrected  V and  V data  obtained  in  this  research. 

P s 

Figures  7 through  9 present  the  variation  of  these  velocities  with  the 
direction  of  wave  speed  measurement , 0.  These  figures  consider  only  the 
tests  at  low  density,  with  figures  for  each  of  the  octahedral  stress 
levels  of  5,  10,  and  20  psi.  Figures  10  through  12  present  the  same 
sequence  of  data,  but  consider  only  the  tests  at  high  relative  density. 

Appendix  B presents  the  same  information  as  in  Figures  7 through 
12,  but  in  the  form  of  tables.  This  appendix  includes  a separate  table 
for  each  test,  presented  in  order  of  test  number.  These  tables  include 
the  composite  correction  factor  resulting  from  the  corrections  described 
in  sections  4.1.1  through  4.1.4,  as  well  as  the  resulting  corrected 
velocity.  These  tables  also  include  the  reference  number  of  the  match- 
ing photo  in  Appendix  A that  documents  the  waveform. 


*t.3  Least-Squares  Curve  Fitting 


I first  attempted  to  fit  straight  lines  and  parabolas  through 
Cartesian  coordinate  plots  of  each  of  the  various  groups  of  data  in 
Figures  7 through  12.  While  both  fit  about  equally  well,  it  soon  be- 
came obvious  that  each  would  produce  a discontinuity  in  V versus  0 when 
0 passed  through  either  0°  or  90°.  This  is  not  compatible  with  the 
physical  situation.  I then  went  to  the  more  suitable  polar  coordinate 
system  and  sought  another  form  of  curve. 

A number  of  engineers  have  used  the  following  form  of  equation, 
adapted  to  the  present  problem,  to  express  anisotropic  behavior  of  soils 

vo  - v90  * (vo  - V °°s29  (2) 


where  numerical  subscripts  represent  degrees.  Casagrande  and  Carrillo 
(19M)  first  intuitively  proposed  an  equation  of  this  type  to  express 
soil  strength  anisotropy.  They  also  noted  that  such  an  equation  ex- 
presses the  strain  in  any  direction  in  terms  of  the  principal  strains  in 
an  elastic  material.  Schmertmann  (196U)  presented  the  same  equation, 
in  a different  form,  to  express  the  induced  variation  in  effective 
stress  on  any  plane  in  an  anisotropic  stress  field.  Lo  (1965)  also  used 
this  form  of  equation  to  express  his  measured  undrained  strength  anisot- 
ropy in  various  clays  and  included  it  in  his  analytical  correction 
procedure  for  anisotropic  strength  effects.  It  seems  well  established 
that  this  form  of  equation  can  be  used  to  express  anisotropic  soil  be- 
havior, both  inherent  due  to  soil  structure  and  induced  due  to  aniso- 
tropic stress.  I therefore  used  Equation  2.  I found  the  standard 
error.  Equation  3,  about  the  same  as  when  I tried  Cartesian  lines  or 
parabolas — thus  further  confirming  the  suitability  of  Equation  2. 

Each  best-fit  curve  is  graphed  in  Figures  7 through  12.  Table  2 
presents  the  results  in  terms  of  the  best-fit  values  of  and  VQ  and 
the  standard  error  of  the  estimate,  S , given  by: 

—I  (3) 
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In  almost  all  cases  we  used  all  the  data  points  available  to  us 
to  determine  the  best-fit  Equation  2.  This  number  of  points  varied  from 
5 to  11.  However,  within  Tests  3 and  10  we  discarded  a few  points  as 
being  nonrepresentative  of  the  remainder  of  the  data. 

The  results  of  the  aforementioned  curve  fitting  also  provide  a 

basis  for  describing  the  degree  of  anisotropic  behavior  of  the  sand 

with  respect  to  V and  V . Consider  the  ratio  for  such  a pur- 

r p s 0 90  * 

pose,  for  both  the  push  and  shear  waves.  Table  3 presents  the  individual 
and  averaged  VQ/Vg0  ratios  as  taken  from  Table  2. 

4.4  Determinations  of  Bulk  Moduli,  B 

One  can  now  use  the  best-fit  curves  illustrated  in  Figures  7 
through  12,  and  summarized  in  Table  2,  to  determine  the  variations  in 
B with  0.  Putting  Equation  2 for  both  and  Vg  into  Equation  1 yields 
an  equation  for  B as  a function  of  0 in  terms  of  VQ  and  V . Use  of  the 
various  VQ  and  values  in  Table  2 then  permits  solving  for  B versus 
0 for  each  density,  K,  and  octahedral  stress  condition. 

For  more  efficient  presentation,  the  writer  chose  to  transform  the 

various  VQ  and  values  at  = 5 and  20  psi  to  equivalent  values  at 

10  psi  and  then  average  these  with  the  real  values  at  10  psi  to  obtain 

an  overall  average  at  10  psi.  For  this  transformation  the  writer 

0.2 

assumed  V aQ‘  , as  described  in  section  4.1.1.  Table  4 lists  these 
overall  average  best-fit  velocities  at  10  psi.  Use  of  the  Table  4 
values  for  and  V^q  in  the  equation  for  B produced  the  results  graphed 
in  Figure  13.  This  figure  presents  one  of  the  final  products  of  this 
research  and  illustrates  the  inherent  sand  structure  and  the  stress 
induced  anisotropy  of  B.  The  values  of  B plotted  in  Figure  13  result 
from  Equation  1 and  involve  the  prior  assumption  of  elastic  isotropy. 

The  resulting  anisotropy  of  B indicates  that  the  assumption  of  isotropy 
is  not  valid  in  our  sand,  especially  when  subjected  to  an  anisotropic 
stress  system. 


5.  DISCUSSION  OF  RESULTS 


5.1  Anisotropy  of  Sand  Structure 

The  heavy  lines  in  Figures  7 through  12  indicate  the  best-fit 
Equation  2 curves  expressing  the  variation  in  either  or  Vg  with  the 
direction  of  wave  velocity  measurement,  0.  These  figures  and  Table  4 
clearly  show  that  these  velocities  do  not  remain  constant  with  direction, 
even  in  the  K = 1 isotropic  stress  field.  There  seems  to  exist  a per- 
haps slight,  but  definite,  trend  for  velocities  to  increase  as  the  wave 
travel  direction  changes  from  the  vertical  (0  = 90°)  to  the  horizontal 
(0  = 0°).  Table  3 shows  VQ/V^0  averages  of  I.067  for  and  1.133  for 
Vg  at  low  density,  and  1.120  and  1.134  at  high  density.  This  probably 
reflects  an  inherent  anisotropy  in  the  sand  structure  itself  at  all 
densities,  resulting  from  its  raining  deposition  in  a gravity  field. 

The  falling  and  then  striking  and  bouncing  sand  grains  probably 
arrange  themselves  with  a statistical  preference  for  their  long  and 
intermediate  axes  to  lie  on  the  horizontal  plane,  which  would  represent 
the  position  of  minimum  potential  energy.  Thus,  waves  propagating  in 
the  horizontal  direction  would  statistically  have  fewer  grain-grain 
gaps/contacts  to  jump  across.  Because  such  jumps  presumably  slow  the 
wave,  it  follows  that  with  such  axis  preference  the  waves  would  travel 
faster  in  the  horizontal  direction.  Hence,  at  K = 1.0  we  found  both 
and  Vg  greater  at  0 = 0°  for  all  y^  and  a\ 

Reference  to  Equation  1 shows  that  the  difference  between  V and 
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V controls  B.  B would  not  necessarily  change  with  0 if  both  V and  V 
s V ° 

change  in  the  same  direct:’  able  3 shows  they  do  so  change.  The 

K = 1 curves  in  Figure  iat  when  considering  the  overall  average 

results  from  this  rei.oa.ot,  one  can  say  that  B is  approximately  iso- 
tropic in  an  isotropic  stress  situation  in  our  pluvially  deposited  sand. 
The  calculations  used  to  plot  Figure  13  show  that  at  low  density 

B0/Bg0  = 1.01  and  at  high  density  BQ/BgQ  = 1.21. 
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5.2  Effect  of  an  Anisotropic  Stress  Condition 


It  seems  apparent  from  Figures  7 through  12  and  Tables  3 and  U 

that  the  shear  wave  velocity  is  not  strongly  influenced  by  an  isotropic 

stress  condition  if  we  make  comparisons  at  the  same  level  of  octahedral 

stress.  The  Table  U average  V at  K = 1 equals  1026  ft/sec  and 

s u 

1059  ft/sec  at  K = 1/3.  The  similar  comparison  for  V ^ is  916  and 
933  ft/sec.  The  differences  fall  within  experimental  error.  An  iso- 
tropic V agrees  with  the  conclusion  of  Hardin  and  Black  (1966)  and  also 
s 

stated  in  Richart,  Hall,  and  Woods  (1970)  that  the  shear  modulus  and 
shear  velocity  at  a point  depend  only  on  the  magnitude  of  octahedral 
effective  stress  and  not  on  the  principal  stress  ratio  at  that  point. 


Contrary  to  the  above,  our  data  show  that  is  measurably  aniso- 
tropic in  an  anisotropic  stress  situation.  The  data  from  the  K = 1/3 

tests  seem  to  show  clearly  that  V decreases  as  0 decreases.  The 

P 

Table  ^ average  V equals  1701  while  equals  lU60  ft/sec,  a dif- 
ference well  outside  experimental  error.  This  means  we  get  the  greatest 
compression  wave  velocity  with  waves  traveling  in  the  same  direction  as 
the  major  principal  stress  (vertical)  and  the  least  velocity  when  travel- 
ing in  the  direction  of  a minor  principal  stress  (horizontal). 


To  further  confirm  the  above,  I also  performed  a supplemental 
Test  13a  in  which  we  reversed  the  direction  of  principal  stresses  and 
performed  a few  check  velocity  measurements  with  K = 2.2.  We  performed 
these  tests  after  releasing  the  stresses  from  Test  13,  which  produced 
some  local  slumping  of  the  sample  around  its  top  perimeter.  We  then  let 
the  sample  sit  over  the  weekend.  Nevertheless,  after  reapplying  aniso- 
tropic stresses  to  bring  the  sample  to  o'q  = 5.0  psi,  but  with  K = 2.2, 
we  now  found  maximum  = 1513  ft/sec  in  the  horizontal  direction  and 

minimum  V = 985  in  the  vertical.  Note  that  this  gives  a V^/V.^  ratio  of 
p u yu 

1.5^  compared  with  an  average  of  0.87  when  K = 1/3.  Figure  7(b)  in- 
cludes the  V „ and  V points  (solid  triangles)  from  Test  13a.  These 
pO  p90 

points  differ  markedly  from  the  K = 1/3  data.  It  seems  well  confirmed 

that  V behaves  anisotropically  in  an  anisotropic  stress  field. 

P 
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The  stress-anisotropic  behavior  of  has  an  exaggerated  effect  on 
B because  of  the  nature  of  Equation  1,  which  involves  the  subtraction  of 
two  large  numbers  to  give  relatively  small  differences.  Figure  13  shows 
clearly  the  stress- induced  anisotropic  behavior  of  B.  B^/B^  has  now 
decreased  to  0.32  at  low  density  and  0.42  at  high  density,  compared  with 
the  corresponding  K = 1 values  of  1.01  and  1.21. 

Note  that  the  ratios  0.32/1.01  = 0.32  and  0.42/1.21  = 0.35  seem 
remarkably  close  to  the  stress  ratio  K = 0.33.  This  suggests  that  B 
depends  in  an  important  way  on  the  effective  stress  acting  on  the  plane 
normal  to  the  direction  of  wave  travel. 

5.3  Dr.  Woods,  Consultant 

As  planned.  Dr.  R.  D.  Woods,  Professor  of  Civil  Engineering  at  the 
University  of  Michigan  and  an  acknowledged  expert  in  the  field  of  soil 
dynamics,  served  as  an  active  consultant  to  this  project.  While  at  the 
University  of  Florida  as  a Fugro  Fellow,  he  helped  with  the  initial 
planning  of  this  project.  He  then  followed  the  progress  of  the  work  by 
means  of  letter  reports  from  the  writer.  He  also  came  to  the  University 
of  Florida  April  14-16  and  spent  two  and  a half  days  working  with  us  on 
the  project.  Before  this  visit  he  reviewed  many  of  our  first-arrival 
interpretations  by  studying  the  then-available  typical  waveform  photos, 
presented  here  in  Appendix  A.  During  this  visit  he  reviewed  in  detail 
all  our  interpretations  for  and  Vs  from  Tests  10,  11,  and  12.  In  all 

this  review  he  found  no  consistent  and  major  differences  between  his  and 
our  interpretations.  He  did  find  minor  differences  and  an  occasional 
major  difference.  But,  even  the  acceptance  of  his  interpretations  would 
not  change  significantly  the  best-fit  equations  given  herein. 

On  April  15,  Dr.  Woods  participated  through  the  entire  Test  13.  He 
did  most  of  the  first-arrival  interpretations.  His  first-arrival  inter- 
pretations resulted  in  velocities  that  fit  in  well  with  our  interpreta- 
tions in  previous  replicate  Tests  3 and  9. 

Dr.  Woods  has  also  reviewed  this  report  and  agrees  in  general  with 


the  results  and  conclusions  presented  herein. 


5.4  Errors  and  Accuracy- 

Table  5 presents  a list  of  all  of  the  known  possible  sources  of 
error  in  this  work,  together  with  some  comments  concerning  the  random- 
ness of  the  error,  its  possible  or  likely  magnitude,  and  the  possible 
importance  of  the  error.  Note  that  in  section  4.1  four  of  these  sources 
of  error  (Nos.  6 through  9)  were  discussed,  and  we  applied  appropriate 

corrections  to  the  raw  data  for  these  sources  before  plotting  V and  V 

P s 

in  Figures  7 through  12. 

By  far  the  most  important  uncorrected  source  of  error  comes  from 
the  difficulty  in  interpreting  the  first  arrivals.  This  proved  only  an 
occasional  problem  with  the  compression  waves,  due  mostly  to  prearrival 
of  the  wave  in  the  second  transducer  as  a result  of  short-circuiting 
around  the  chamber  steel-water-rubber  perimeter.  More  serious  was  the 
more  common  problem  cf  detecting  the  first  arrival  of  the  shear  wave. 
Appendix  A presents  a typical  set  of  waveform  data  for  each  or  Vs 
determination  and  the  reader  can  judge  for  himself  or  herself  about  the 
interpretation.  As  noted  in  section  5-3,  Dr.  Woods  agreed  with  the 
interpretations  he  checked. 

Note  that  our  values  are  too  high  compared  with  predictions  from 
the  well-known  empirical  equations  from  Hardin  and  Drnevich  (1972),  also 
found  in  Richart,  Hall,  and  Woods  (1970),  p.  154),  as  well  as  the  more 
recent  ones  from  Kuribayashi  et  al.  (1975*  unpublished).  Our  shear 
velocities  are  approximately  30%  higher  than  predicted  using  the  Hardin- 
Black  equation  and  about  20 % higher  when  using  the  newer  Kuribayashi 
equation.  This  may  indicate  an  error  in  our  results  or  that  the 
equations  do  not  apply  to  our  sand  and/or  our  chamber  test.  Both  of 
these  equations  result  from  interpretations  of  laboratory  data  using 
only  resonant  column  testing. 

The  writer  looked  into  the  possibility  that  different  magnitudes 
of  shear  strain  might  account  for  the  above  differences  in  shear 
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velocity.  Using  equation  5 for  shear  strain,  y,  suggested  by  Dr.  Woods 
and  the  second  form  of  it  that  results  from  assuming  simple  harmonic 
motion, 


Y 


A 

■V  -PP  - i\j 

- V T - 
s s 


(5) 


where 

A = maximum  peak  to  peak  amplitude 
Vs  = velocity  of  the  shear  wave,  ft/sec 

T = period  of  shear  wave  = l/f  sec 

S 2 
a = maximum  acceleration  ft/sec 
max  2 

g = gravitational  constant  = 32.2  ft/sec 
f = frequency,  Hz 

__o 

the  writer  obtained  a shear  strain  range  of  about  0.05  to  0.5  x 10  % at 

the  first  arrival  of  the  shear  wave  for  a representative  group  of  10 
determinations  during  this  research.  These  strain  magnitudes  fall  below 

_3 

the  10  % noted  by  Hardin  and  Drnevich  as  below  those  for  maximum  shear 

modulus,  Gmaxs  and  therefore  fall  in  the  range  for  which  they  report 
their  equation  applicable.  Dr.  Woods  reported  the  data  from  Kuribayashi 
et  al.  involved  strains  of  10  S».  Differing  strain  levels  do  not 
appear  to  explain  the  differences  in  V . 

In  addition  to  the  shear  wave  velocities,  the  ratio  of  V /V  mea- 

P s 

sured  in  this  research  in  the  K = 1 tests  varied  from  about  l.U  to  1.6, 
suggesting  a Poisson's  ratio  of  about  0.0  to  0.2.  This  seems  too  low 
compared  with  the  values  of  0.25  to  0.U  usually  suggested  for  sand. 
However,  the  0.0-0. 2 range  results  from  computations  assuming  an  iso- 
tropic, elastic  half-space,  which  we  found  did  not  exist  even  with  the 
K = 1 stress  condition. 

5.5  Supplemental  Test  1^  for  Vg 

Because  of  a nagging  uncertainty  in  the  accuracy  of  our  Vg  inter- 
pretations, we  decided  to  perform  a check  test  using  a more  accurate 

2k 


method  for  detecting  the  first  arrival  of  the  shear  wave.  Adapting  the 
work  of  Anderson  and  Woods  (1975).  we  made  a miniature  scissors-type 
wave  generator  which  we  then  used  to  produce  inverted  polarized  shear 
waves  (Figure  14).  This  device  expands  out  and  slides  down  when  struck 
downwards,  and  again  expands  out  but  slides  up  when  struck  upwards.  This 
produces  the  same  compression  wave  but  a different  polarization  of  the 
shear  wave  when  alternately  struck  up  and  down.  The  point  of  divergence 
in  overlaid  successive  waveforms  gives  the  point  of  shear  wave  arrival. 

We  used  this  technique  successfully  in  supplemental  Test  ll+. 

Test  lU  repeated  the  previous  Tests  3,  9,  and  13  with  the  same  dry 
sand  at  the  low  density  condition  and  the  boundary  stresses  at  K = 1. 
However,  because  we  had  sent  two  accelerometers  away  for  service,  we  had 
only  two  available  for  Test  1**.  We  placed  these  at  the  A1  and  A2  loca- 
tions, giving  0 = 30°,  and  introduced  the  waves  at  the  side  access  hole 
(see  Figure  6).  After  some  experimentation  with  the  scissors  device, 
we  were  able  to  get  reproducible  divergent  waveforms  to  clearly  define 
the  first  arrival  of  the  shear  wave.  The  final  photos  in  Appendix  A 
(Nos.  218  through  220)  document  the  typical  sets  of  superposed  waveforms 
we  obtained  in  supplemental  Test  l4,  and  our  interpretations  for  Vg 
(and  V ) first-arrival  times.  Appendix  B also  includes  a table  for  the 

correction  factors  and  the  resulting  final  average  V and  V values  we 

s p 

obtained  from  Test  l4. 

The  open  triangles  in  Figures  7(a),  8(a),  and  9(a)  show  the 

Test  lk  results  for  Vg  together  with  the  previous  values  from  our  other 

methods  of  interpretation.  The  solid  triangles  in  these  figures  show 

our  Test  ll+  values  for  V . The  reader  will  see  that  these  additional 

P 

values  for  Vg  fit  well  with  our  previous  determinations  at  the  10-  and 
20-psi  stress  levels  and  appear  superior  to  our  previous  determinations 
at  0 = 30°  at  5 psi,  which  we  had  already  discounted  as  anomalous  and 
likely  in  error.  Test  l4  further  confirmed  the  overall  validity  of  our 
previous  V interpretations  wherein  we  did  not  have  the  advantage  of 
inverted  shear  waveforms. 
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6.  CONCLUSIONS 

6.1  - The  scope  and  accuracy  of  the  present  research  appears  adequate 

to  answer  the  research  question  concerning  the  possible  anisotropy 
of  the  bulk  modulus  in  an  isotropic  stress  field. 

6.2  - The  compression  and  shear  wave  velocities  vary  with  direction  in 

an  isotropic  stress  field  due  to  inherent  anisotropy  in  the  sand 
structure.  While  observable,  these  trends  seem  to  have  only  a 
modest  (0-20%)  effect  on  producing  an  inherent  anisotropy  in  B. 

6.3  - This  work  tends  to  confirm,  within  +10  percent,  that  the  shear 

wave  velocity  depends  only  on  the  level  of  octahedral  effective 
stress  and  not  on  stress  anisotropy. 

6.4  - In  contrast  with  V , the  compression  wave  velocity,  V , increases 

s p 

significantly  in  the  direction  of  the  major  principal  stress. 
depends  primarily  on  the  effective  stress  on  the  plane  normal  to 
wave  travel. 

6.5  - The  values  of  B determined,  listed,  and  plotted  in  this  report 

result  from  the  use  of  Equation  1 and  involve  the  prior  assumption 
of  elastic  isotropy.  The  resulting  anisotropy  of  B shows  that  this 
assumption  and,  therefore.  Equation  1 are  not  correct.  A correct 
equation  would  require  the  use  of  at  least  two  values  each  of  an 
elastic  modulus  and  Poisson's  ratio  to  express  a minimum  degree  of 
elastic  anisotropy. 
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TABLE  1 - CHRONOLOGY  OF  TEST  OPERATIONS  FOR  A TYPICAL 
TEST  3 TO  13 


Day 

Operation 

Time 

Required 

Men 

Required 

1 

Empty  chamber  from  previous  test, 
weigh  sand  as  removed,  put  into 
dryer  if  needed 

3 hr 

2 

Set  up  chamber  for  next  pluvial 
deposition 

3 

2 

2 

Rain  sand  to  form  specimen  and 
place  four  accelerometers  at 
proper  positions 

3 (low) 

Yd 

7 (high) 

Y , 

1 or  2 

Seal  chamber  and  complete  test 
setup 

'd 

3 

1 

3 

Apply  K = 1 or  K = 1/3  chamber 
pressures  to  5-psi  test  level 

iA 

1 

Hold  pressure  for  15  min 

lA 

Alternately  get  V and  V from 
each  pair  of  accelerometers , 
do  five  times , get  a typical 
scope  waveform  photo 

1-1/2 

Repeat  above  for  other  two  pres- 
sure levels 

It 

Release  pressures,  close 
down  chamber 

1 

2 

I 


Average  total  hours  = 21 
Average  total  man-hour  = h5 


TABLE  2 - RESULTS  OF  THE  LEAST-SQUARES  FITTING  OF 
EQUATION  2 THROUGH  CORRECTED  DATA 
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V 

P 

ve 
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(vo  - V 
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cos  6 
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or 

No. 

V90 

v0 
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(Psi) 
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vg 

Pts 

(ft/sec 

) (ft/sec 

) (ft/sec) 

(ft/sec 

1 Fig. 

5 

1.0 

V 

p 

11 

1106 

1269 

1.1U7 

62 

7a 

V 

s 

6 

806 

578 

0.717 

40 

1/3 

V 

p 

6 

1422 

1272 

0.895 

67 

Tb 

V 

s 

5 

888 

978 

1.101 

37 

10 

1.0 

V 

p 

11 

1325 

1375 

1.038 

111 

8a 

V 

s 

9 

853 

975 

1.143 

43 

1/3 

V 

p 

6 

1619 

1369 

0.846 

35 

8b 

V 

s 

6 

878 

ll44 

1.303 

28 

20 

1.0 

V 

p 

11 

1513 

1538 

1.017 

127 

9a 

V 

s 

9 

1000 

1122 

1.122 

67 

1/3 

V 

p 

6 

1791+ 

1556 

0.867 

71 

9b 

V 

s 

.6 

978 

1069 

1.093 

27 

5 

1.0 

V 

p 

7 

1253 

11*38 

1.148 

50 

10a 

V 

s 

6 

825 

866 

1.050 

42 

1/3 

V 

p 

9 

1606 

1328 

0.827 

97 

10b 

V 

s 

7 

897 

928 

1.035 

103 

10 

1.0 

V 

p 

7 

1453 

1609 

1.107 

50 

11a 

Vs 

6 

9U7 

1106 

1.168 

57 

1/3 

V 

p 

9 

1791 

1503 

0.839 

52 

lib 

V 

s 

7 

997 

1000 

1.003 

114 

20 

1.0 

V 

7 

1666 

18U1 

1.105 

60 

P 

12a 

V 

s 

6 

1091 

1291 

1.183 

44 

1/3 

V 

9 

2013 

1775 

0.882 

62 

P 

V 

s 

7 

1116 

1353 

1.212 

139 

12b 

TABLE  5 - POSSIBLE  SOURCES  OF  ERROR  IN  V AND  B DETERMINATIONS 


No. 

Type  Error 

Estimated 
Maximum 
% Effect  on 

V V 

_J3  _s_  Both 

Random 

Comments 

1 

First  arrival  interpretation 

5 25 

No 

Some  independent 

checks  by  Dr.  Woods 

2 

Accelerometer  vibrations  not 
fixed  in  surrounding  sand 

Unknown, 
but  prob- 
ably small 

No 

Accelerometer  4 cm 
long,  2 cm  test 

wave-lengths  5 to 

15  cm,  c-c  spacing 
35  to  69  cm 

3 

Estimated  times  on  scope 

4 

Yes 

Reduces  to  2%  when 
average  5 

1* 

Inaccurate  accelerometer 
placement 

2 

Yes 

Placed  to  +1/4  cm 

5 

Nonuniform  stress  condition  at 
top  chamber  due  to  rigid  top 
platen 

Unknown 

No 

Affects  only  0=0° 
and  30°  tests,  tend 
to  increase  Vp  at 

K = 1/3,  tend  to 
reduce  anisotropy 
of  B 

6 

Varying  levels  of  average  a' 
between  accelerometers 

3 

No 

Corrected  (see  4.1) 

7 

Scope  sweep  rate 

0.5 

No 

Calibrated  and 
corrected 

8 

Sand,  compression  under  chamber 
stress 

1 

No 

Measured  and 
corrected 

9 

Varying  average  whole-chamber  y 
between  tests  in  a density 
group 

d 1 

Yes 

Measured  and 
corrected 

10 

Relatively  loose/dense  layers 
between  accelerometers 

1/2 

Yes 

Minor  due  to  repli- 
cate tests  and 
averaging 

11 

Loosening  dense  sand  and 
densifying  loose  sand 
around  accelerometer 
during  placement 

1/4 

No 

12 

Different  levels  of  shear 

strain  during  V and  V 
# P s 

determinations 

Probably 

small 

No 

Approximately  50$  of 
push  impact  energy 
used  in  shear  im- 

pacts,  both  strain 
less  than  10-3? 


FIGURE  1 - SCHEMATIC  OF  UNIVERSITY  OF  FLORIDA  CALIBRATION  CHAMBER  AND  LOADING  SYSTEM 
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(c)  W.  Whitehead  positioning  A3 


FIGURE  3 - RESEARCH  EQUIPMENT 
AND  METHODS 


d)  H.  Margas  (L)  and  Dr.  Woods 
(R)  interpreting  waveforms  during 
Test  13 
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(a)  Photo  shows  the  vertical 
rods , attached  to  the 
3-in.  diameter,  rough- 
bottomed  plates  buried 
3 cm  below  sand  surface,  passing  through  the  two  holes  in  the  top  platen. 
They  are  encased  in  an  insulating  plastic  pipe  to  prevent  shorting  of  the 
signal  trigger  system.  The  cables  from  the  four  buried  accelerometers 
come  out  through  the  larger  center  hole.  All  accelerometer  cables  were 
carried  as  far  as  possible  away  from  wave  paths  between  accelerometers 
within  the  sand  sample 


(b)  Ball  pulled  back  to  reference 
bar,  then  let  go  to  swing  and 
strike  horizontal  bar  clamped 
to  vertical  rod.  Weight  was 
added  to  insure  no  plate  slip 
during  impact 


(c)  Ball  dropped  from  reference 
bar  height  to  strike  center 
of  plate 


FIGURE  4 - THE  MECHANICAL  SWINGING  AND  DROPPING  BALL 


IMPACT  SYSTEM  USED  IN  TESTS  7 THROUGH  13a 


I 

£ 

FIGURE  5 - ACCELEROMETER  PATTERN  IN  PLANE  THROUGH  AXIS  OF  CHAMBER,  USED 

FOR  TESTS  3 AND  U 


INITIAL  BOTTOM  SAND  I WATER  BAG  LOAD  CELL  ON  TOP  PISTON 

t 


FIGURE  6 - ACCELEROMETER  PATTERN  IN  PLANE  THROUGH  AXIS  OF  CHAMBER,  USED 

FOR  TESTS  6 THROUGH  13 


0 <100  800  1200  1600  2000  ft/sec 

FIGURE  9 - BEST-FIT  EQUATION  2 CURVES  THROUGH  V_  AND  V_ 

P s 

DATA  AT  LOW  DENSITY,  OCTAHEDRAL  STRESS  =20.0  PSI 


(a)  K = 1 

tests  6,11 


800  1200  1600  ft/sec 


(b)  K = 1/3 
tests  7,10,12 


800  1200  1600  ft/sec 


FIGURE  10  - BEST-FIT  EQUATION  2 CURVES  THROUGH  Vf 

AND  V DATA  AT  HIGH  DENSITY,  OCTAHEDRAL  STRESS 
s 

= 5.0  PSI 
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Notes:  a)  Used  Fig.  5 accelerometer  pattern 

b)  Sand  wet  after  test  due  to  membrane  leakage.  Unii 
wts.  estimated  from  prev.  tests,  prob.  high. 


TEST  NO. 


0?  = 5.0  psi  j 
iTd  = 92.3  pcf  I 
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Notes:  a)  Used  Fig.  5 accelerometer  pattern 
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Notes:  a)  First  test  with  Fig.  6 accelerometer  pattern. 
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TEST  NO.  ? 


Notes:  a)  Minor  membrane  leakage  found  after  test.  Wt.  dry 
sand  corrected  by  estimating  wt.  water  (30  lb), 
b)  H.  Margas  did  1st  arrival  interpretations. 
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Notes:  a)  First  test  to  obtain  V at  9 = 0°  by  striking  side 
of  tank  with  ball. 
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TEST  NO.  !0 


K = 'A 
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Notes:  a)  V data  involving  A4  (9  = 90°,  60°)  suspect  because 
does  not  fit  previous  data  patterns.  Later  careful 
removal  of  A4  found  no  obvious  reason  for  problem, 
b)  Velocity  interpretations  checked  by  Dr.  Woods. 
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Notes:  a)  Velocity  interpretations  checked  by  Dr.  Woods. 


B5 


TEST  NO.  U K = '/) 


i 

ii  - 

• i 

o 

psi 

o 

- 10.0  psi 

%'= 

20.0  psi 

tfd 

= tots 

pcf 

i* d = /Ol.i  pcf 

<Td  =/02.t  pcf 

wave  4 
with 
horiz. 

correct 

factor 

corr. 

1^. 

)hotc 

No. 

correct 

factor 

corr. 

vp  vs 

foto 

No. 

correct 

factor 

corr. 

foto 

No. 

0 = 

.399 

IW 

It) 

E 

If? 

.3f3 

i»n 

Ifi 

90° 

.9  IT 

f)  3 

m 

■9 

nr 

.393 

1291 

iro 

.993 

Df) 

nr 

Mi 

mo 

i?i 

.593 

it  ft 

10/ 

60° 

1.001 

rgi 

its 

.359 

sn 

nr 

.5^ 

5 If 

102 

1.0)3 

/2)s 

ni 

1.013 

mo 

m 

/.003 

lf)l 

IfV 

30° 

1.0)3 

1)0 

/tt 

1.015 

if) 

t.o  or 

uoo 

if) 

n 

/. 01 ? 

/2SI 

no 

I.OOi 

/ns 

in 

.333 

l?Ol 

106 

Notes:  a)  Velocity  interpretations  checked  by  Dr.  Woods. 
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Peak  to  peak  maximum  vibration  amplitude 
Maximum  acceleration  in  simple  harmonic  motion 
Dynamic  hulk  modulus  or  modulus  of  volume  compressibility 
B when  0=0° 

B when  0 = 90° 

Frequency,  Hz 
Maximum  shear  modulus 

2 

Gravitational  constant  = 32.2  ft/sec 
Principal  stress  ratio  = a^/a^ 

Standard  error  (see  Equation  3) 

Period  of  the  shear  wave 

Velocity  of  the  compression  or  push  wave 

Velocity  of  the  shear  wave 

Velocity  of  push  or  shear  wave  when  traveling  in  0 direction 
VQ  when  0=0° 

VQ  when  0 = 90° 

Shear  strain 
Dry  unit  weight 

Angle  between  direction  of  wave  travel  and  plane  on  which 
acts 

Mass  density  = Y^/g 

Octahedral  normal  stress  = (a^  + + a^)/3 

Major  principal  normal  stress 

Intermediate  principal  normal  stress,  = this  research 
Minor  principal  normal  stress 


NOTE:  All  primed  qualities  used  for  stresses  are  understood  to  be 
effective  stresses. 
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